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Abstract. - The wavefunction of an incommensurate ground state for a one-dimensional dis- 
crete sine-Gordon model - the Frenkel-Kontorova (FK) model - at zero temperature is calculated 

eby the quantum Monte Carlo method. It is found that the ground state wavefunction crosses 
over from an extended state to a localized state when the coupling constant exceeds a certain 
' ' critical value. So, although the quantum fluctuation has smeared out the breaking of analyticity 

' , transition as observed in the classical case, the remnant of this transition is still discernible in 

H , the quantum regime. 

o 



The Frenkel-Kontorova (FK) model was proposed more than half a century ago to study 
commensurate-incommensurate phase transitions . It has since become a widely used model 
not only in condensed matter physics but also in nonlinear dynamics (see e.g. In condensed 
matter physics, it serves as a model for an adsorbate layer on the surface of a crystal j^, Q , ionic 
conductors and glassy materials , charge-density waves [|| , and chains of coupled Josephson 
junctions 0. Moreover, the FK model is also a nonlinear lattice model having a normal 
thermal conductionJ|g] . Recently there has been a growing interests in applying this model to 
study dry- friction ||^|and atomic scale friction or nanotribology In nonlinear dynamics, 
the FK model provides a vivid physical realization of the seminal Kolmogorov-Arnol'd-Moser 
theorem. 

The FK model describes a chain of particles connected by harmonic springs and subject 
to an external sinusoidal potential. The classical behavior of this model has been studied 
by Aubry |^. As a discrete sine-Gordon model, the FK model exhibits many new features 
that are absent in the continuum sine-Gordon model. One of the most striking features is 
the so-called "breaking of analyticity" transition. In an incommensurate ground state, as the 
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coupling constant increases beyond a critical value, the FK model displays a transition between 
the unpinned state (sliding phase) and the pinned state (pinned phase). This transition is 
characterized by the appearance of a phonon gap, the discontinuity of the Hull function, and 
etc. These phenomena have been also found in many generalized FK models pH |. 

Although a lot of studies have been done on the classical FK model, very little study has 
been devoted to the quantum FK model ^ Such study is also of more than 

academic interest. For example, recently there has been a surge in the study of atomic scale 
friction, or nanotribology. This problem is of tremendous importance to technology. The FK 
model has been used to study this problem. However, in the nano-regime, it's inevitable that 
quantum mechanics will play an important role. So an understanding of the quantum FK 
model will be important to the understanding of nanotribology. 

In this Letter we would like to study the quantum version of the one dimensional FK model 
by using the quantum Monte Carlo approach. We shall concentrate on the wavefunction 
of an incommensurate ground state. As we shall see later the wavefunction changes from an 
extended state to a localized one as the coupling constant is increased beyond a certain critical 
value. This crossover is a quantum remnant of the classical breaking of analyticity transition. 

The Hamiltonian operator of the one dimensional FK model is 



2m 2 



(1) 



Here m is the mass of the particles, 7 the elastic constant of the spring, 27r/go the period of the 
external potential, V the magnitude of the external potential, and a the equilibrium distance 
between two nearest neighbor particles. For convenience, we shall rescale the variables so as 
to obtain a dimensionless Hamiltonian 



(2) 

L 

where K = Vq^/^ is the couphng constant. We define an "effective Planck's constant": 

(3) 




This effective Planck's constant is the ratio of the natural quantum energy scale (?iwo) to the 
natural classical energy scale (7/90) (iHj ^^^^ = 7/m. 

To study the ground state of Hamiltonian (H), the quantum Monte Carlo method will be em- 
ployed [p^ . As in the classical case, we will concentrate our attention on the incommensurate 
state characterized by the golden mean value era = (V5 — 1)/2. In the classical case the winding 
number is defined as the average separation of atoms per period, i.e. a = limAr_>oo '^2Ntt" • 
a is a rational number, the ground state is commensurate, and if it is an irrational number the 
ground state is called incommensurate. In the quantum case, Xi{i = 0, 1, • • • , N) is defined 
as the expectation values of the particles. As usual, we use the method of continued fraction 
expansion and approximate ug by its rational convergents uq = Fn/ Fn+i, where Fn is the nth 
Fibonacci number defined by the recursion relation, Fn+i — Fn + Fn^i with _Fo = 0,i^i = 1. 
Therefore, in our quantum Monte Carlo computation, we choose Fn+i particles which are 
embedded into Fn external potentials with period of 27r. We impose periodic boundary 
condition on the chain. 

Since the external potential is periodic and has period 27r, we can fold the wavefunction 
to this period and then take the average over all particles in the interval [0, 27r]. This 
quantity gives the probability of finding particles at a given potential position X. We plot 
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the averaged probability density (l^'P) in Fig. 1 with 144 particles embedded in 89 potentials 
for a fixed h (= 0.2) for different values of K. (Notice that (|vl'P) is normalized, namely, 
/o = 1.) We observe that, in the smaU K regime, the probability of finding the 

particles at any place of the potential is almost the same (see the curve for K = 0.1.). This is 
quite similar to the sliding phase in the classical version. However, as the coupling constant 
increases, the probability of finding the particle at the top of the potential decreases, while 
that at the lower part of the potential increases, as is noticeable from the appearance of peaks 
in the curves. As the potential goes beyond a certain critical value, the probability of finding 
the particles at the top is almost zero, see e.g. the curves ioi K — 2 and 5. 

In the classical FK model. Coppersmith and Fisher |l^ have proposed a "disorder parame- 
ter" to describe the transition from the pinned phase to unpinned phase. This parameter is de- 
fined as the minimum distance of a particle from the top of a well, Del = Hiiuj-.n \Xj''—2n{n+^)\. 
It can be seen that, as long as Dd ^ 0, the particles are pinned, while Dd — corresponds 
to an unpinned state. This "disorder parameter" also measures the discontinuity (or width 
of the biggest gap) of the Hull function. One might naively try to use the same function 
to describe the quantum crossover. For instance, one may define a very similar quantum 
disorder parameter Dq = minj_„ \Xj — 27r(n -I- i)|, where is the expectation value of the 
j/'th particle's position. However, this parameter Dq could not capture the crossover of the 
ground state wavefunction. The reason is that, in the quantum case, the particle can tunnel 
from one side of the potential to the other, thus the gap in the classical Hull function does not 
survive the quantum fiuctuations (see Refs. and [|l^). It turns out that even if K > K^, 
Dq might be also close to (or equal to) zero. Therefore, a new parameter is needed to describe 
the wavefunction crossover in the quantum FK model. To this end, we define the probability 
of finding particles at the potential top as a quantum "disorder parameter" , which is denoted 
as Ff. Pt is given by 

N 



1 / 1 

P*=nT. \^iX)\'s(X-2n{n+- 



2) ) dX. (4) 

In Fig. 2(a)-(b) we plot Pt versus K for different temperature and different system size. 
In Fig. 2(a), we fix the quantum fiuctuation 'h{= 0.2), and the particle numbers (13/21) by 
changing the temperature from Tg — 0.2/30,0.2/120, to 0.2/480. The convergence is quite 
fast as is seen from the figure. The sharp decrease of Pt is very evident for K between 1 and 
2. In this small K regime, Pt changes dramatically: it drops almost two orders of magnitude. 
This dramatical change can also be seen from Fig. 2 (b). There the temperature is fixed at 
Tf. = fi/NSr = 0.2/120, which can be regarded as effectively zero temperature. The system 
size is changed from 21 particles to 144 particles. In the quantum Monte Carlo method the 
Feynman path integral presentation is used. In this presentation, solving a one dimensional 
quantum problem is equivalent to solving a two dimensional effective classical problem with 
an extra dimension as the imaginary time JTtI . The size of this extra dimension is fi(3^ where 
/? — 1/Tf,, Tf, is the temperature. Therefore, in the zero temperature limit, we shall take [3 
infinite. 

It is worth pointing out that another parameter which can be used to depict this crossover 
is the maximal fiuctuations of the particles. In the quantum case, we have observed taht the 
particle situated at a position nearest to the maxim of the potential has always maximal 
fluctuation since it has the largest classical potential energy. This observation has been 
verified numerically by the quantum Monte Carlo results, and theoretically by our squeezed 



state theory |15 . Thus, we can use this maximal fiuctuation as another measure of quantum 



"disorder" . To make a qualitative comparison with the classical disorder parameter Dd defined 
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Fig. 1. - The wavefunction of an incommensurate ground state having winding number 89/144 at 
fixed h = 0.2 for different values of K. The wavefunction (l^'j^) means the probability of finding the 
particles at X (mod 27r). The curves are for K — 0.1,1,2, and 5, respectively. The wavefunction 
becomes localized at the lower part of the potential as K changes to 2 and 5. 

by Coppersmith and Fisher ||l^, we shall make use of the standard deviation, i.e. the square 
root of the maximal fluctuation, 

A = max{((X2)-(X,)2)^}. (5) 

Here the average (• • •) is taken over all the paths produced in the quantum Monte Carlo 
simulation, which is about 4,000 in our calculations. A versus K is plotted in Fig. 3. The 
computations given in this figure have been carried out with cr = Fn/Fn+i = 34/55. 

It is noticeable that the transition of the wavefunction is characterized by the different 
i^-dependent behavior of A. For K < K*, A is a constant. It does not depend on K, but 
it changes with fi. For K > K*, the maximal fluctuation increases with K, but it does not 
change with h. Furthermore, this quantum parameter A is approximately equal to the classical 
disorder parameter Dd- For comparison with the classical disorder parameter, we also include 
Del in the inset of Fig. 3. In the classical case, the transition occurs a.t K = K*, where 
K* = 0.971635.... 

These results show that the width and the shape of the probability density do not change by 
changing the external potential strength in the small K(< K*) regime. The width only depends 
on the strength of the quantum fluctuation h. However, this picture changes dramatically when 
K > K*. In this regime, the proflle of the probability density spreads out, and the width of 
the probability density is insensitive to quantum fluctuations; instead it depends only on the 
coupling constant. It must be stressed that the analogy between Dd and A cannot carried 
too far, since A is not exactly the quantum correspondence of Dd- This is why even if we let 
h go to zero, A would not be zero in the regime oi K < K*. 

In summary, we have studied the ground state wavefunction of the FK model at zero 
temperature. The wavefunction undergoes a crossover from an extended state (analogous to 
the sliding phase) to a localized state (analogous to the pinned phase) as the external potential 
increases. Therefore, although the quantum fluctuation has smeared out the breaking of 
analyticity transition as observed in the classical case, the remnant of this transition is still 
discernible in the quantum regime. 

We would like to thank S. Aubry, F. Borgonovi, D. K. CampbeU, H. Chen, R. B. Grifhths, 
H.-Q. Lin, L.-H. Tang, and W.-M. Zhang for helpful discussions during several periods of the 
work. BL thanks the Abdus Salam International Centre for Theoretical Physics (Trieste) for 
its kind hospitality during his visit there in the summer of 1998. This work was supported in 
part by the grants from the Hong Kong Research Grants Council (RGC) and the Hong Kong 
Baptist University Faculty Research Grant (FRG). 
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Fig. 2. - The quantum "disorder parameter" Pt, i.e. the probability of finding particles on the top 
of the external potential versus K for a fixed quantum fiuctuation h = 0.2. (a) Pt for different 
temperature at a fixed winding number 13/21. The temperature changes from 0.2/30,0.2/120, and 
0.2/480, respectively, (b) Pt for different winding numbers: 13/21, 34/55, and 89/144, at fixed 
temperature T = 0.2/120, which can be effectively regarded as zero. 



Fig. 3. - A versus K. Different symbols represent different effective Planck's constant, which are 
?i = 0.01, and 0.1, respectively. We draw the classical disorder parameter defined by Coppersmith and 
Fisher [17] in the inset for comparison. The winding number for all calculations given in this figure is 
34/55. 
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